Huge (but finite) time scales in slow relaxations: beyond simple aging 
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Experiments performed in the last years demonstrated slow relaxations and aging in the con- 
ductance of a large variety of materials. Here, we present experimental and theoretical results for 
conductance relaxation and aging for the case-study example of porous silicon. The relaxations are 
experimentally observed even at room temperature over timescales of hours, and when a strong 
electric field is applied for a time the ensuing relaxation depends on t^. We derive a theoretical 
curve and show that all experimental data collapse onto it with a single timescale as a fitting param- 
eter. This timescale is found to be of the order of thousands of seconds at room temperature. The 
generic theory suggested is not fine-tuned to porous silicon, and thus we believe the results should 
be universal, and the presented method should be applicable for many other systems manifesting 
memory and other glassy effects. 

PACS numbers: 71.23.Cq, 72.15.Cz, 72.20.Ec, 71.23.-k 



For a large variety of systems in nature, the relaxation 
back to equilibrium depends strongly on the system's his- 
tory. Examples of such behavior are abundant in the 
magnetic properties of spin-glasses [ll-Q , the mechanical 
properties of polymers and biological material 0-0] , and 
in the behavior of colloids near the jamming transition 
• Understanding the aging properties is important for 
applied physics and engineering, as well as for theoretical 
physics, as this phenomenon is interrelated to the physics 
of glasses. 

In the last few decades much progress has been made 
in understanding aging behavior in these broad range of 
physical systems, both theoretically and experimentally 
[8|. In one of the common aging protocols, the system is 
perturbed out of equilibrium for a time tu,, the waiting- 
time. A time t after the perturbation is switched off, a 
physical observable in the system is measured. It was 
found that for various systems, the signal depends only 
on the ratio t/t^, which is referred to as 'full' or 'simple' 
aging. Examples of full aging exist for a diversity of phys- 
ical systems, such as electron glasses @-12|, spin-glasses 
[Tst . vortices in superconductors and granular sys- 
tems Nontheless, full aging is not generic. In some 
cases, a heuristic form of scaling is used fjt/tlf;,) [lij . 
which was recently given theoretical support [l7[. Many 
systems, however, show more complex behavior, which is 
the case discussed here. 

In this work, we study, experimentally and theoret- 
ically, the aging properties of a system of porous sil- 
icon, whose conductance exhibits aging behavior over 
timescales ranging from seconds to days. It should be 
emphasized that the theoretical model does not con- 
tain ingredients which are specific to porous silicon, and 
should be universally applicable to a broad range of 
complex systems whose dynamics is governed by ther- 
mal activation or equivalent mechanisms. While some 
of the measurements appear to be close to full aging. 



clear and significant systematic deviations from it arise 
when the involved experimental timescales increase. This 
does not conform to the /{t/f^) scaling form mentioned 
above. The results suggest that the system must have 
unique, system-dependent finite timescale, characteriz- 
ing the breakdown of the full aging regime. Here, we 
consider a model that predicts a particular form of aging, 
and demonstrate how it fully accounts for the complex 
form of aging in porous silicon. We succeed in accurately 
deducing a timescale characterizing the system, and find 
that, experimentally, it is sensitive to the temperature as 
well as to the time since the sample has been manufac- 
tured, i.e., its 'age' in the usual sense. 

The structure of the manuscript is as follows. We 
first present the porous silicon system and the experi- 
mental protocol. We proceed by outlining the physical 
model, which was used by some of us earlier in the con- 
text of electron glasses, and put forth the prediction for 
the aging behavior, including the explicit dependence on 
the system dependent timescale. This timescale has the 
physical significance of the reciprocal of the lower cut- 
off of the underlying relaxation distribution. Finally, we 
show the excellent agreement of experimental data and 
theory, both for the 'excitation' and 'relaxation' parts of 
the experimental protocol, which will be presently de- 
fined. 

Experimental System.- There has been much interest 
in the physics of porous silicon, due to its potential in 
applications such as biosensors and its remarkable 
optical properties [l^, and due to the high controlla- 
bility of its fabrication. Mesoporous silicon (MesoPS) 
consist of a disordered network of interconnected silicon 
nanocrystals. For the samples studied here the size of 
the nanocrystals is of the order of tens of nanometers 
[lit . The samples are fabricated by electrochemical etch- 
ing of heavily doped p-type silicon in hydrofluoric acid 
solution. The linear dimensions of the samples are of 
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the order of tens of microns. Since the obtained mesoPS 
layers are not thin, the role of contact resistance and 
of the silicon substrate in the electrical characterization 
is negligible. The electronic conduction mechanism in 
this material is not fully understood. Previous studies 
have suggested that the conduction is limited by various 
narrow channels of porous silicon, where due to the pres- 
ence of nearby trapped charges the electronic passage is 
severely suppressed However, the conduction de- 

pendence on temperature in some cases is reminiscent of 
variable-range- hopping [2l| . which may point out to ex- 
istence of other possible conduction mechanisms at lower 
temperatures. Full details of the experimental system are 
presented in |22|. 

The I-V characteristics of the porous silicon samples 
are empirically found to be superohmic. This implies 
that when pushed out-of-cquilibrium by an electric field, 
the time dependent conductance {i.e., the current to volt- 
age ratio) must increase to its higher out-of-equilibrium 
value. Remarkably, it does so very slowly, over the course 
of hours. 

The protocol used in the experiment is as follows. 
First, the sample is let to relax for a long time, of the 
order of several hours, in a constant electric field F, pre- 
sumably reaching very close to its steady-state. The 
value of the field is then raised, for a time t^, which 
can be anything from seconds to a few hours. The ini- 
tial value of the electric field is now restored, and the 
system relaxes to its initial steady-state. The current is 
measured as a function of time, through the experiment. 
Fig. [T] shows the results for a particular value of t^. 

Model- In it was argued that various, generic 



physical mechanisms may lead to a broad distribution 
of relaxation rates which is approximately: -P(A) ^ 1/A. 
One of these, which we believe to be the relevant one 
here, relies on thermally activated processes. Indeed, if 
we assume such processes, A ^ e^'^/^, and a distribu- 
tion of barriers which is approximately uniform between 
a lower and upper cutoff, P{U) « C, we obtain the above 
distribution of relaxation rates (also between a lower and 
upper cutoff, Amin and Amax)- This simple physical rea- 
soning is at the heart of seminal works on 1 // noise [23| . 

In the following we will adopt the aging picture de- 
scribed in Ref. [2J], with an important difference: here, 
the cutoffs of the distribution will play a major role, 
and wc will not assume that we are in the asymptotic 
regime where they do not govern the physics, namely 
l/Amaa; <C t <C 1 / Xmin wiU uot bc valid hcrc. Let us 
review the basic mechanism: when the external pertur- 
bation (voltage in our case) is applied, the slow modes 
are partially excited, and the conductance of the system 
increases to a larger out-of-equilibrium saturation value. 
In complete analogy to charging a capacitor, the exci- 
tation of a single mode with characteristic frequency A 
is [1 — e^'^*]. Using the above distribution for the slow 
modes, we thus obtain: 



(5cr^^(£) 



Amin 



-[l-e-^*]dA, 



(1) 



with i the time elapsed from the application of the per- 
turbation. 

Thus in the excitation regime: 

,5^'=" (i) ^A + Ei [A,nax<] - E, [A„,i„il , (2) 
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FIG. 1: Measurement of time-dependent conductance. The 
electric field exerted on the system, F, is changed from an 
initial value to a larger value for a time t^, as shown in the 
figure. Upon the increase of field, the conductance increases 
monotonically, attempting to reach a new, higher steady-state 
value, but does not saturate on the experimental timescales. 
Upon the restoration of the electric field to its initial value, 
the conductance relaxes in a form different from exponential. 



with A = log(Amax/Amin)! SLild Ei thc exponential in- 
tegral function. As long as we do not consider too short 
times, i.e., as long as i ^ 1/Amaa;, we can neglect the 
term Ei [Amaxi] • 

Similarly, in the aging regime (see Fig. [T|), each of the 
modes relaxes exponentially with rate A. Therefore the 
excess conductance is: 



(5cr"s(t,i™) (X / -[l-e-^*"']e-^*dA. (3) 

Which can be approximated for times which are longer 
than 1/Ainax as: 



5a''^{t,t^)(xEi[\^,^t]-Ei[K,,,,{t + t^)]. (4) 

Were we to assume that the experimental timescales 
are much shorter than the reciprocal cutoff Amin, we could 
approximate the exponential integral function by a log- 
arithm. Under that approximation, the cutoff drops out 
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of the final expression for the excess conductance, yield- 
ing the form log(l + t^/t). This is the case for electron 
glasses, which has been experimentally observed for a 
variety of systems The remarkable thing about the 
porous silicon system, is that the timescale associated 
with the cutoff Amin is long enough such that one can 
measure slow relaxations over the course of hours, but, 
it is still experimentally measurable, and its finite value 
must be considered in order to understand the a-priori 
complex aging behavior. In fact, the deviations from full 
aging in this case arise precisely from the finite value of 
Amin (and will disappear as the experimental timescales 
will become shorter and shorter). 

Results. - Fig. [2]compares Eq. ^ and the experimental 
results. Amin is taken as a fitting parameter. The excel- 
lent agreement using a single fitting parameter to data 
collapse five curves, shows that the theoretical framework 
suggested here is adequate. It thus allows us to associate 
a single timescale T^ax = 1/Amin with each porous sili- 
con sample (which is typically of the order of thousands 
of seconds), which has the physical meaning of being the 
longest relaxation time existing in the system. 

Fig. [3] shows the analysis of the same data set in an- 
other way, by rescaling time according to Eq. Q , leading 
to a linear dependence of the excess conductance on the 



3 

D 

D 2 













i„ = 6000s 




— i„ = 3000s 




— i„ = 1200s 




— i„ = 500s 




— i„ = 300s 


\ \ \^ 


Theory with 

T,nax = 4800s 




---log[l + t^/t] 







10' 



t/tnj 

FIG. 2: Experimental results and theory compared, for a spe- 
cific porous silicon sample, measured at T = 299K. The x 
axis shows the time divided by the waiting-time tn,, which 
is different for each of the curves. The y axis shows the 
excess conductance from its equilibrium value. The ex- 
periment was performed for waiting-time values of t^n — 
300s, 500s, 1200s, 3000s, 6000s. Full aging corresponds to the 
case were there would be complete data collapse of all curves. 
Here, while the curves for the smaller values of are close 
to complete data collapse, the curves corresponding to the 
larger waiting-times are significantly off from it. The solid 
lines correspond to the theoretical form of Eq. Q, with a 
single fitting parameter Amin, for all of the five curves. 



rcscaled variable, i.e., we present a parametric plot where 
the y-axis is the measured conductance and the x-axis is 
given by Eq. (g]). 

Having determined the timescale Tmax from the ag- 
ing experiments, we test its applicability also during the 
excitation step, where the conductance is expected to in- 
crease as an exponential integral function, as described 
by Eq. ([2]). Fig. ^ shows the good agreement between 
experiments and theory, where now Amin is no longer 
a fitting parameter. This is consistent with the sug- 
gested framework, in which the modes which are excited 
during the application of the electric field are the same 
ones which relax during the aging phase. 

So far we have established the physical significance 
of a timescale characterizing the system. Having done 
that, one can now investigate the dependence of Tmax on 
the system's parameters {e.g: temperature), and learn 
about the underlying physical mechanisms. For the case 
of porous silicon, we have found that Tmax is sensitive 
to the temperature, which is consistent with the ther- 
mal activation mechanism. For the sample measured, 
upon changing the temperature from BOOK to 280K the 
timescale Tmax changed from about 800 seconds to more 
than 5000 seconds (it should be noted that the measure- 
ments of Figs. [2][4]wcre done on a different sample, giving 
different timescales) . Using measurements at several dif- 
ferent temperatures would allow one to find the cutoff 
of the barrier distribution, i.e., the deepest traps avail- 
able. Moreover, we have found that leaving the samples 
exposed to air for a long period (of the order of several 
weeks) substantially inffuences this timescale: the value 
of Tmax was 160 scconds for a freshly made sample, but 
grew to 4800s for a two months old sample. This is an- 



4 

^3 

3 

D 2 



1 



1 2 





= 6000s 




= 3000s 




= 1200s 




= 500s 




= 300s 


3 


4 



El [Xmin t] — El [Xmin {t + tyj)] 

FIG. 3: The same data of Fig. [2]is analyzed in a different way; 
the time (x axis) is rescaled according to Eq. leading to 
a linear dep endence of the excess conductance. For electron 
glasses [24l [. rescaling according to log(l -I- tm/t) would be 
adequate, and one would not be able to use this procedure to 
deduce a timescale. Here, the exponential integral function is 
necessary, and not only its asymptotic form. 
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other form of 'aging', which is maybe closer to the every- 
day use of this terminology: as the sample grows older, 
its typical timescale grows. Here, it means, remarkably, 
that the effective traps in the porous silicon are not static 
in time, but are 'deepened' as time progresses. This is 
a further process, occurring on much longer timescales 
than those involved in the experimental protocol, which 
justifies the analysis with a constant in time T,-aax- We be- 
lieve our analysis demonstrates the power in the method, 
which allows one to gain insights into the microscopies of 
porous silicon, through the macroscopic measurements of 
the conductance's slow relaxations. For example, the de- 
termination of the system's characteristic timescale was 
recently discussed in the context of electron glasses [lEl , 
and the method presented here might prove useful also 
in that case. 

Summary.- We studied slow relaxations and aging in 
a system of porous silicon, which is well suited for this 
purpose since it supports relaxation modes which are ex- 
tremely long, yet it has a finite, measurable timescale 
Tmax- We have shown that one can understand quan- 
titatively the aging behavior of the conductance in the 
system, and data collapse the measurements from five 
different experiments using a single (sample dependent) 
fitting parameter, the timescale T„iax- The underlying 
theory used was not particular to porous silicon, and 
certain limits of it, where the aging is independent of 
the timescale r^nax, have been shown to describe other 
systems [131 1 suggesting that the results are robust. We 
have shown that r^ax is temperature sensitive, consistent 




FIG. 4: The change in conductance during the 'excitation' 
phase, when the system is driven out-of-equilibrium, for the 
same sample and conditions of Figs. [2] and [31 In accordance 
to Eq. ((2]), the rise in conductance is described by an ex- 
ponential integral function, with a single fitting parameter 
Amin. Its value was determined by the aging characteristics 
(see Fig. [2]), and therefore in the above figure Amin is not a 
fitting parameter. The good agreement between the theory 
and experiment prove that the excitation of the conductance 
results from the same physical origin as the relaxation. 



with the picture of thermal activation over large energetic 
barriers in the systems. Since the theoretical model does 
not contain any ingredients particular to porous silicon, 
we believe these intriguing phenomena should be appli- 
cable to various other systems where deviations from full 
aging are observed, in material science and in other glassy 
systems. 
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